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The paper concerns Weyl-Heisenberg covariant SIC-POVMs (symmetric infor- 
mationally complete positive operator valued measures) and full sets of MUBs 
(mutually unbiased bases) in prime dimension. When represented as vectors 
in generalized Bloch space a SIC-POVM forms a, (P — 1 dimensional regular 
simplex (d being the Hilbert space dimension). By contrast, the generalized 
Bloch vectors representing a full set of MUBs form d + 1 mutually orthogonal 
d — I dimensional regular simplices. In this paper we show that, in the Weyl- 
Heisenberg case, there are some simple geometrical relationships between the 
single SIC-POVM simplex and the d -I- 1 MUB simplices. We go on to give 
geometrical interpretations of the minimum uncertainty states introduced by 
Wootters and Sussman, and by Appleby, Dang and Fuchs, and of the fiduciality 
condition given by Appleby, Dang and Fuchs. 
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1. Introduction 
There has been much interest in recent years in SIC-POVMs P, ^ ^, 

1221 1331 131] (symmetric informationally complete positive operator valued measures; 
citations in order of first appearance online or in print). SIC-POVMs have been 
constructed analytically in Hilbert space dimension d — 2-13, 15 and 19 (existence 
of analytic solutions for d = 11, 15 communicated to author privately |35j : for rf = 15 
also see ref. [36 ), and numerically in dimension 5-45. The fact that they exist in 
every dimension up to 45 (at least to a very high degree of numerical accuracy) 
means one may plausibly speculate that they exist in every dimension — although 
this remains to be proved. 

MUBs [ini [in [H Ell Eg [Ml HOI Hlllia 113111113 (mutually unbiased 

bases; only a representative selection of papers cited) have been the subject of 
intense investigation for a rather longer period of time. It is known that the number 
of MUBs in dimension d cannot exceed d+1, and that the maximum number of d+1 
MUBs exist whenever c? is a power of a prime number. However it remains an open 
question whether the maximal number of d-l- 1 MUBs exist in any non-prime-power 
dimension. 

The purpose of this paper is to describe some geometrical relationships between 
these somewhat elusive structures when the dimension d is a prime number, so 
that full sets of d-l- 1 MUBs certainly exist, and SIC-POVMs exist at least in some 
cases. We also give a geometrical interpretation of the minimum uncertainty states 
introduced by Wootters and Sussman |13 HI] and Appleby, Dang and Fuchs jSU] , 
and of the related fiduciality conditions given by Appleby, Dang and Fuchs [30] . 
Our analysis builds on previous work by Bengtsson and Ericcson [101 112j and Ap- 
pleby [17]. 

We win work in generalized Bloch space [Tni[ni[13[5ni[nTl[51[S3[Sl[53[n3[57]. 

If d = 2 the Bloch representation of an arbitrary density matrix p is defined by 

P^lil + h-a) (1) 

where b is a 3-vector having length < 1 (the Bloch vector) , and cti , a2, era are the 
Pauli matrices. Pure states correspond to the case |b| = 1, and mixed states to 
the case |b| < 1. The crucial feature here is that b • cr is a trace-zero matrix. So if 
d > 2 we can proceed in an analogous manner, by writing 

P=^(l + B) (2) 

where B is a trace-zero Hermitian matrix acting on c?-dimensional Hilbert Space Hd- 
More technically, B is an element of the Lie algebra s\i{d) (although we will make 
no use of the theory of Lie algebras in this paper). To obtain the analogue of the 
Bloch vector b some authors proceed by choosing a basis for su(d) (for example, the 
Gell-Mann matrices). We, however, find it more convenient to regard the matrix 
B itself as a vector, and to define an inner-product on su(d) by the formula 

(B,B') = -^tr(BB') (3) 

(readers familiar with the theory of Lie algebras will perceive that, apart from a 
scale factor, this is just the Killing form on su((i)). Note that this gives us a real 
inner product space. We define 

||B|| = v/(B:B) (4) 



to be the corresponding norm. 
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If d = 2 the set of Bloch vectors corresponding to density matrices is very easily 
described: it is just the ball ||B|| < 1. If d > 2 the geometry is rather more 
intricate. Let B be the Bloch body (the set of Bloch vectors corresponding to 
density matrices). Then [50l|53l[56j B is wholly contained by the out-ball 

6o = {B e su(d): ||B|| < 1} (5) 

and wholly contains the in-ball 

6, = {Besu(d): IIBII <l/(d-l)} (6) 

More succinctly: 

BiCBCBo (7) 

Moreover Bo and Bi are, respectively, the smallest and largest balls with these 
properties. In other words 

(1) For each < r < 1 there exists a Bloch vector B € S such that ||B|| = r 

(2) For every r > l/{d ~ 1) there exists a vector B such that ||B|| = r and 
B(jtB 

The set of pure states is the intersection of B with the out-sphere So = {B G 
su(ci):||B||=l}. 

In terms of this picture an orthonormal basis in the Hilbert space Hd corresponds 
to a d— 1 dimensional regular simplex in Bloch space, wholly contained in the Bloch 
body and with all of its vertices lying on the out-sphere So- To be more specific, 
let \tpi), . . . , \ipd) be an orthonormal basis in Hd and let B^ — d\ipa){'>Pa \ — 1 be the 
corresponding Bloch vectors. Then 

d 

J2^a^0 (8) 

a=l 

and 

(Ba,Bb) = ^^(dSab-l) (9) 

Now let . . . , be a second orthonormal basis, and let B'j^, . . . ,BJ; be the 
corresponding Bloch vectors. The two bases are mutally unbiased with respect to 
one another if and only if | {ipal^'b) I ~ ^ ^'^^ ^- terms of Bloch space this 
means the bases are mutually unbiased if and only if the corresponding simplices 
are mutually orthogonal: 

(B.,B',)=0 (10) 

for all a,b. It follows that a full set of d + 1 MUBs [TOl [H HZ] corresponds to a 
family of d -|- 1 mutually orthogonal d — 1 dimensional simplices in Bloch space. 
They thus define a highly symmetric polytope in Bloch space. This picture also 
makes it clear |T2 [57] why the question, whether full sets of d + 1 MUBs exist 
in every dimension d, is hard to answer. The problem is that the out-sphere So has 
dimension d^ — 2, whereas the manifold of pure states BDSq has dimension 2d — 2, 
which is smaller if d > 2, and very much smaller if d is large. It is easy to construct 
a family of d -f 1 mutually orthogonal d — 1 dimensional simplices whose vertices lie 
on the out-sphere So, but if d > 2 it is much less easy then to rotate this structure 
so that its vertices all lie on the lower dimensional sub- manifold B O So- 

A SIC-POVM [Tni[Ill[57], by contrast, corresponds to a single d^ — 1 dimensional 
regular simplex in Bloch space, wholly contained in the Bloch body, and with its 
vertices all on the out-sphere So- To be more specific, let |V'i), . . . IV'd^) be the 
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vectors defining a SIC-POVM. So 

1 

^Y.\^r){A\=0 (11) 

and 

Ki^rmf + (12) 

d + 1 

for all r, s. In terms of the corresponding Bloch vectors this means 

^B^ = (13) 

r=0 

and 

(B„B,) = ^^(d2^,,-l) (14) 

for all r, s. It is hard to prove that this structure exists in dimensions > 2 for essen- 
tially the same reason that it is hard to prove full sets of MUBs exist in dimensions 
> 2: namely [lOl [El [27j , though it is easy to construct a — 1 dimensional regular 
simplex with its vertices all on the out-sphere So, it is very difficult then to rotate 
the simplex so that its vertices all lie on the lower dimensional sub-manifold BOSo- 
The purpose of this paper is to examine the relationship between the (P — 1 
dimensional simplex corresponding to a SIC-POVM, and the family of d-l-l mutually 
orthogonal d ~ 1 dimensional simplices corresponding to a full set of MUBs. We 
will confine our attention to the case where d is a prime number > 2 (the case 
d = 2 being almost trivial). We will further restrict our attention to SIC-POVMs 
and MUBs which are covariant under the action of the Weyl-Heisenberg group. 

2. The Weyl-Heisenberg and Discrete Symplectic Groups 

We begin by enumerating some salient facts concerning the Weyl-Heisenberg and 
Discrete Symplectic groups which will be needed in the sequel. For the proofs of 
the statements which follow the reader may consult, for example, refs. [TSllSO]. 

To define the Weyl-Heisenberg group choose an orthonormal basis |0), . . .\d— 1) 
for Hd and define operators X and Z by 

X\a) = \a+1) (15) 
Z|a)=cj"|a) (16) 

for all a, where addition of indices is modulo d and where w = g27ri/d^ rpj^^ Weyl- 
Heisenberg displacement operators are then defined by 

where the integers qi, q2 range over the interval [0, d — 1] and where denotes the 
unique integer in the range [0, d — 1] such that 2^^ x 2 = 1 mod d (so, for example, 
= 3 if d = 5). Note that these definitions depend on the fact that d is odd; for 
even d the definition of the displacement operators is a little more complicated. 

We will also have occasion to use the discrete symplectic group SL(2,Zrf), con- 
sisting of all matrices of the form 

^=(; ?) 
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where a, (3, 7, 5 are integers modulo d such that det F = 1 mod d (and where we are 
again making use of some simphfications due to the fact that d is odd). For given 
F define a corresponding unitary operator Up by [15] 

(where, as before, the notation (2/3)"^, 2~^ signifies the inverse in the modular 
sense; it should also be noted that we are making use of the fact that d is specifically 
an odd prime — if d is non-prime the definition of Up is a little more complicated). 
Then 

UpD^Ul = Dp^ (20) 

for all p. 

We can obtain a full set of MUBs by acting on the standard basis |0), • • • ,\d—l) 
with symplectic unitaries Up- For the sake of definiteness define 

G=(j ;) (21) 



and 



I Uh\o) m — 00 

Then it is readily verified that the bases |m, 0), |m, 1), . . . |m, rf — 1) as to ranges 
over the set {0, 1, ... (d — 1), 00} constitute a full set of d + 1 MUBs. Note that the 
basis for which to = is just the standard basis. The choice of the label to = 00 
rather than to = d to denote the last basis is motivated by the discrete affine 
plane picture . It is shown in ref. [301 that this is the only set of MUBs which 
can be obtained by acting on the standard basis with symplectic unitaries Up (up 
to permutation and re-phasing). To describe the action of the Weyl-Heisenberg 
displacement operators it is convenient first to define operators X^n, Z„i by 

X„|TO,a) |TO,a+ 1) (24) 

and 

Z„i\m, a) — uj'^ "'\m,a) (25) 
(c./ Eqs. ([15|) and p6|l ). In terms of these operators we then have 

D„\m, a) = < „-i , , (26) 

(in other words the Weyl-Heisenberg displacement operators act on each basis by 
permuting and rephasing). 

A SIC-POVM covariant under the action of the Weyl-Heisenberg group is defined 
by first choosing a fiducial vector 1-0) with the property 

mD^m' = ^{d5^^o + l) (27) 
for all q (assuming this to be possible), and then defining 

IV'q) = D^W (28) 
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3. Geometrical Relationships of the SIC and MUB Simplices 

These preliminaries concluded we are now in a position to analyze the geometrical 
relationships between a single SIC simplex and the d + 1 MUB simplices (in odd 
prime dimension, for Weyl-Heisenberg covariant SICs and MUBs). 

Let B™ ,j be the Bloch vector corresponding to \m,a){m,a\, and let be the 
Bloch vector corresponding to \'4'q){ipq\- So [TUllT^[?fj 

(B-.,B-,,J J^('^^--l) (29) 

and 

(B^,B^,) = ^^(dX.q'-l) (30) 

Let Vm denote the subspace spanned by the vectors B™ q, . . . B™ It will also 
be convenient to choose an orthonormal basis for Bloch space. Let B™ ^ „ and „ 
be the components of B^ ^ and B^ respectively relative to this basis (where the 
index u runs from 1 to — 1). 

Now define Q„j to be the matrix with elements 



^m,a,u^m,a,u' (31) 



a 

a=0 

It is easily seen that Qm is the projection operator projecting onto the m*'^ MUB 
hyperplane Vm- 

Q™B-,. = jf- (32) 
ID m m 

We wish to calculate the projections of the Bloch vectors B^. We have 



,B 



•1 d 

a=0 



-^tr(B-,ByB-„ (33) 



d2 

a=0 

Now 

tr(B^°_,B^) = tr((rf|77i,a)(m,a| - l) (dji^qXV'ql - 1 

— d^l (to, a|'(/'q)|^ — d 
= d^ I (to, aj-DqlT/il^ — d 

_ I d'^P7n,(a-qi+mq2) - d TO 7^ OO 

\d^Pni,{a+q2)~'d TO = OO 

where we used Eq. (|26p in the last line, and where 

Pm,a = \ {m,a\%l^)\^ (35) 

is the probability of obtaining outcome a when measuring the TO*'' MUB in the 
fiducial state \'ip). Taking into account the fact that 



(34) 



Y.^Z^a = Q (36) 



a=0 
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for all m we deduce 



Q^W^={^\+'-'''' ' (37) 



m = oo 



where 



C„i,a — \/<j + 1 Pni.,f)B„i,b+a (38) 



f)=0 

It was shown in ref. [301 that 



Pm,bPm,b+a = -7— — (^a.O + 1) (39) 

. a+1 



b=0 



Consequently 



(C,„^a, Cm,a') = idda,a' - 1) (40) 

from which it can be seen that the vectors Cm, a constitute a d— 1 dimensional regu- 
lar simplex in the MUB hyperplane Vm- This gives us two regular simpliccs in each 
MUB hyperplane: the one comprising the vectors B™^, and the one comprising 
the vectors C„^. 

Putting all this together we conclude 

(1) The vector makes the same angle with each of the hyperplanes Vm- 
To put it another way: the Bloch vectors representing the SIC-POVM are 
maximally distant from the MUB hyperplanes. 

(2) For each to, the vectors Bq constituting the SIC-simplex project onto 
only d vectors. More specifically, if m 7^ oo the d vectors {B^^^^^^ ^^y. b = 

0, 1, . . . , c? — 1} project onto the single vector {d + l)~^/^Cm,a, while if 
m = 00 the d vectors {B^j^ 5 = 0, 1, . . . , d — 1} project onto the single 

vector {d+l)-^/'^C„„a- 

(3) For each to the d — 1 dimensional projection of the d^ — 1 dimensional SIC 
simplex is itself a regular simplex. 

Since the vectors B™ ^ and Cm.a both constitute regular simpliccs of the same 
size it must be possible to rotate one onto the other. We now investigate this 
rotation. 

Let Rm be the matrix with elements 

^ _ -[^ rf-i 

Rm,uv = ^ ^ Cm.a.uBjj^ g^ y (41) 

It is easily seen that 

RmRm = Qm (42) 

and 

^mB™ ,j = C„ (43) 

So the restriction of Rm to Vm is a real orthogonal transformation taking the B™ ^ 
simplex onto the a simplex. 

We next diagonalize Rm, regarded as a complex matrix. Define 

BZa^^^^J2^'""^^,b (44) 

6=0 
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for a = 1, . . . , d — 1 (the definition does, of course, also make sense for o = 0; 
however B™ o = 0)- It is straightforward to verify 

B;5i,„ = ^=L=Ea;-«^B-, (45) 
(where a runs from to d — 1) and 

d-l 

E *<M = '5a,6 (46) 

«=0 

(where a, b run from 1 to d — 1). So the vectors ^ are an orthonormal basis for 
Vmi regarded as a complex vector space. They are also eigenvectors of Rm'- 

,3 d-l 



\ , ,ctc/^ r>m r>m 

■^m m,a y ^2 / ^ ^m,b,u m,b,v m,c,v 



6,c,u=0 
■ d-l 



,u 



6=0 

d-l 



^ ^ Pm,cBm,c+b,-i 



6,c=0 

-1 d-l 



6,c=0 e=l 
— Pm,aJ->m,,a,u 

where 

d-l 

= Vd+l E ^~"'^Pm,b (47) 

6=0 

Regarded as a complex matrix is unitary, so we can write 

Pm,a = e*^-'" (48) 
for suitable phase angles 9m,a- The fact that „ = Pm,-a means we can assume 

^m^—a ~ 9m,a (49) 

for all m, a. The matrix elements of Rm can then be written 

d-l 

a=l 

We can now use this expression to write Rm explicitly as a rotation matrix in 
canonical form. Define 

= ^ {Bm,-a + B„,) = y^^^ E COS (^) B-, (51) 
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for a = 1, . . . , (d — l)/2. By construction the vectors B^j^,B^j^ are a real or- 
thonormal basis for Bloch space: 

(B^,Q, f,) = (BJ„ BJ„ b) = (5a,b (53) 

(Bf„^„B5„,,) =0 (54) 

for a, 6 = 1, . . . , (d — l)/2. Moreover 



Rm,uv — y ] f cos 0„j_a ^i?„j_a,u^m,a,u "I" ^m,a,u^m,a,v^ 



So when it is expressed in terms of the basis B^ Bj„ . . . , B^ (d-i)/2' ^In (d-i)/2 
the operator i?,„ is the direct sum of the 2x2 rotation matrices 



^cos^m.a -sin6',„,Q 
I sin 6*^, a cos0„i^Q 



(56) 



It foUows from aU this that if we knew the angles 9m,a then we could easily 
reconstruct the SIC-POVM. All we would need to do is to use the operators Rm to 
rotate the MUB simplices B™ ^ onto the projected SIC simplices C„ We could 
then calculate the SIC-POVM using the formula 

= ( Coo,-g2 + X/ ^ni,(qi-mq2) j (57) 

When we embarked on this investigation we entertained the (faint) hope that 
the angles 6m. a would turn out to have a simple form (for instance, that they would 
all be rational multiples of tt). Unfortunately this does not seem to be the case. 
So the result does not seem to take us any closer to solving the existence problem. 
Nevertheless, it does, perhaps, have some intrinsic interest. 

4. Minimum Uncertainty States and the Fiduciality Condition 

The results in the last section provide some insight into the geometrical signif- 
icance of the minimum uncertainty condition introduced by Sussman and Woot- 
ters [IHl |49] and Appleby, Dang and Fuchs ^30] . It was shown by these authors that 
a pure state {ip) minimizes the quadratic Renyi entropy averaged over a full set of 
MUBs if and only if 

a=Q 

for all m, where Pm,a = ajV')]^. 

Let B be the Bloch vector corresponding to an arbitrary normalized state 
Then 

d-l 

QmB=J2Pm,aB^^la (59) 



a=0 



Consequently 




\Q„,B\\' = ^{y^pl^]--^ (60) 



from which it follows that 



(61) 
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for all m if and only if 

d-i 2 

T.Plr,a = J-l (62) 

for all TO. In other words is a minimum uncertainty state if and only if its Bloch 
vector makes the same angle with all the MUB hyperplanes. 

Appleby, Dang and Fuchs 30J further show that is a fiducial state for a 
Weyl-Heisenberg covariant SIC-POVM if and only if it is a minimum uncertainty 
state and, in addition, 

d-l ^ 

^PmMPrn,a+b = -J— J (63) 
a=0 

for all TO and all 6 7^ mod d. This statement also has a simple geometrical 
interpretation, as we now show. 

For each to define Sm to be the operator whose matrix elements are 

Sm,uv — y ] ^m,a+l,u-^m,a,v (64) 

It is readily confirmed that 

^m^m ~ Qm (65) 

So the restriction of Sm to the to*'' MUB hyperplanc Vm is an orthogonal operator. 
It is in fact a rotation operator (as follows from the fact that — Qm which, 
since d is odd, means det(S'„i) = 1). Its significance is that it effects the action of 
the Weyl-Heisenberg displacement operators on Vm- 

DqQmBDl^ = i (66) 

[Sm QmB TO = OO 

Now let B be the Bloch vector corresponding to a minimum uncertainty state. We 
have 

, d-l d-l 

QmB = ^^(B-,„B)B-,, = ^p^UB^.a (67) 

a=0 a=0 



, d-l 

a 



and, consequently, 

(Q„B, S:,Q„,B) = ^ E PM^.b+a - (68) 

fc=0 

So eq. is satisfied if and only if 

(g„B, S:,Q^B) = {ddafl - 1) (69) 

(i.e. if and only if the vectors S','^Q„jB form a d — 1 dimensional regular simplex 
for each fixed to). 

To summarise: a state jV') with Bloch vector B is a fiducial state for a Weyl- 
Heisenberg SIC-POVM if and only if 

(1) B makes the same angle with each MUB hyperplane 

(2) For each to the rotated vectors Q„jB, S'^^Q^^B, . . . , Sf~^Q„^B form a. d—1 
dimensional regular simplex. 
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